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We calculate the energy spectrum of a k = 2 flux tube that is closed around a spatial 
torus, as a function of its length /. We do so for SU(4) and SU(5) gauge theories in 2 



space dimensions. We find that to a very good approximation the eigenstates belong to the 
irreducible representations of the SU(iV) group rather than just to its center, Zn. We obtain 
convincing evidence that the low-lying states are, for I not too small, very close to those 
^ ■ of the Nambu-Goto free string theory (in flat space-time). The correction terms appear to 
be typically of 0(1) in appropriate units, much as one would expect if the bosonic string 
model were an effective string theory for the dynamics of these flux tubes. This is in marked 
contrast to the case of fundamental flux tubes where such corrections have been found to 
be unnaturally small. Moreover we find that these corrections appear to be particularly 
small when the 'phonons' along the string have the same momentum, and large when their 
momentum is opposite. This provides information about the detailed nature of the interactions 
in the effective string theory. We have searched for, but not found, extra states that would 
arise from the excitation of the massive modes presumably associated with the non-trivial 
structure of the flux tube. 



1 Introduction 



Consider aD = 2 + l SU(iV) gauge theory on a three-torus, and a loop of confining flux that 
winds once around a spatial torus of length I. If this flux is in the fundamental representation 
then one discovers [U, [2] , that the low- lying energy spectrum of such a loop is the same as 
that of a corresponding closed string in the Nambu-Goto free string theory (in flat D = 2 + 1 
space-time) up to very small corrections, and that this is so even when I is not much larger 
than the critical length at which the flux loop dissolves. This type of calculation is intended 
to provide detailed information about the effective string theory describing the gauge theory 
and should be useful for theoretical approaches from the string side. 

In the usual field-theoretic picture of the confining string as some non-Abelian version of 
the flux tube in a (dual) superconductor, with some complicated internal spatial structure 
which should provide a complex set of internal modes to excite, it is not easy to understand 
how a closed flux tube that is not much longer than it is wide, can display the phonon-like 
excitations of a thin string to such a good accuracy. It is therefore tempting to see this as 
evidence for some underlying duality where the starting point in describing a flux tube is 
always a thin string, irrespective of its length, just as one has in the usual gauge/gravity 
duality. (See [3] for some recent reviews.) Even in the latter framework, however, one does 
expect states in the spectrum that are additional to those of the free bosonic string and which 
arise, for example, from the non-trivial background metric near the 'horizon' which leads to 
linear confinement. Finding such extra states would provide a useful clue about the nature 
of the hypothetical dual theory. So far no such 'extra' states have been identified in our 
lattice calculations [2] but this may be because the basis of operators we use in our variational 
calculation has little overlap onto any such states (as discussed below). This is one reason 
why it is useful to explore the spectrum of the quite different flux tubes that we study in this 
paper. 

In this paper we explore what happens to the spectrum of a winding flux tube when the 
flux is in a higher representation than the fundamental. In particular for N > 3 there exist 
new, completely stable flux tubes called /c-strings. A fc-string is a flux tube that emanates 
from a source which transforms as ip — > z k ip under a global gauge transformation belonging 
to the center of the group: z e Z^. One may think of this source as being a composite of 
k fundamental sources, together with any number of adjoint sources which do not feel the 
centre. (It is clear that k and k' = N — k are the 'same' and only k < N/2 is interesting.) 
Since k is unchanged by screening with gluons, the lightest flux tube in this sector will be 
stable. That is what one usually calls a fe-string. (Sometimes we shall refer more loosely to all 
such flux tubes, stable and unstable, as fc-strings.) While such a stable flux tube is of especial 
interest, excitations thereof, including flux tubes in other representations, will become stable 
at larger iV and are also of interest. 

Given that a single fundamental flux loop is so well described by a free bosonic string 
theory, one might have expected that a fc-string would simply look like some linear combination 
of free fundamental (k = 1) strings. So, for example, a k = 2 string might look like a 
combination of a k = 1 string that winds twice around the torus, and a state consisting of 
two k — 1 strings, each winding once around the torus. The free string theory would, in this 
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way, predict a simple relation between the string tensions: = koj. However it has been 
known for some time [I] that the lightest /c-string is in fact strongly bound (for k/N not too 
small) and only acquires the free-string value at N = oo (as it must, at fixed k). Given that 
the binding is strong, e.g. in SU(4) one finds au=2 — 1.34a - / <C 2cr/, one can ask to what 
extent such a flux tube also behaves like a free bosonic string and what the corrections to 
such a behaviour tell us about the dynamics. This binding cannot be readily incorporated 
into the usual analytic frameworks [9j [12] since it presumably involves the exchange of small 
contractible loops of string (glueballs) that are not under analytic control. Thus such flux 
tubes are of particular interest. 

The ground state energy of a closed k = 2 flux tube of length I has been calculated in 
(5], where it was found to be accurately described by the free Nambu-Goto string, with string 
tension a = (7fc = 2, and with corrections (in inverse powers of 1/cr/ 2 ) that have coefficients 
that are 0(1) in natural units. This is what one would expect if the k = 2 flux tube were 
to have an effective bosonic string description at large distances, and it provides a contrast 
with the fundamental flux tube where the corrections to the free string description are very 
much smaller PQ[2]. In practice, however, the ground state provides only a limited insight into 
the nature of the string spectrum and so, just as for the fundamental string, it is important 
to perform accurate calculations of at least a few excited states. This is difficult in a lattice 
calculation, where one can only hope to obtain accurate energies for those states which have a 
very high overlap onto the basis that one is using. For the fundamental flux loop we attacked 
this problem, with apparent success, by using a very large basis of operators [2] . Here we shall 
do the same for the k = 2 spectrum, using the basis described in Section [2] and in Appendix 
A. Since we have seen from earlier calculations of this type [H [2], [5j [6] that lattice spacing 
corrections are small, we shall not perform a continuum extrapolation, but will simply choose 
to perform our calculations at fixed, small values of the lattice spacing a. We shall perform 
calculations in SU(4) and SU(5) at (3 = 50 and (3 = 80 respectively. 

In the next section we briefly describe the standard aspects of our lattice calculation and 
then, in some more detail, our choice of operators for winding flux tubes. We follow that up 
with a summary of the corresponding free string Nambu-Goto spectrum, since that will play 
an important part in the interpretation of our results. We then describe our results for the 
spectrum and find that we can interpret at least some of the states within the Nambu-Goto 
picture, simply on the basis of the way the energy depends on the length. We then turn to 
the question of whether the flux tubes know about the full SU(iV) group or just about the Z N 
center and demonstrate that in fact it is the former. We then return to the very interesting 
question of whether the excited states without an obvious Nambu-Goto interpretation might 
not be revealing the excitation of new, internal degrees of freedom. To address this question 
we develop a heuristic method for comparing the shapes of the k = 2 states with the k = 1 
states. This study leads us to infer that all the k = 2 energy eigenstates over which we do 
have control, are most likely to be Nambu-Goto states with, in some cases, large corrections 
to the energy dependence. Comparing those states for which such corrections are small, with 
those for which they are large, enables us to say something qualitative and striking about 
where these corrections arise in the effective string theory. We end with a summary and our 
conclusions. 
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2 Methodology 



Our lattice framework is standard. We shall therefore be brief, except with those aspects that 
are important for understanding the statistical and systematic errors of the present calculation. 

2.1 lattice field theory 

We work on periodic cubic L x x L y X L t lattices with lattice spacing a. The degrees of freedom 
are SU(iV) matrices, Ui, assigned to the links I of the lattice. Our partition function is 



where U p is the ordered product of matrices around the boundary of the elementary square 
(plaquette) labelled by p. The continuum limit is approached by tuning (3 — > oo and in that 
limit (3 = 2N/ag 2 where g 2 is the usual continuum coupling, which has dimensions of mass 
in D = 2 + 1. Our numerical calculations use a standard heat bath/over-relaxation Monte 
Carlo, where one updates all the SU(2) subgroups of the SU(iV) matrices. 

2.2 calculating energies 

We typically wish to calculate the energy spectrum of some sector of states with particular 
quantum numbers (spin, momentum, parity, ...). We do so from correlation functions of opera- 
tors 4>(t) with those same quantum numbers. In a gauge-invariant calculation, the elementary 
components of such operators are typically traces of closed loops, where contractible loops 
typically project onto glueballs and non-contractible (winding) loops typically project onto 
winding loops of flux. In general we have 



C(t) = (0 f (t)0(O)) = ^|(n|0(O)|mc)|V^ ^ |(O|0(O)Hc>|V^ ; E t < E l+l (2) 



where on a lattice we write t = an t , so we see that all the energies will be obtained in lattice 
units, as aEi (n t is the separation in the Euclidean time direction in lattice units). Clearly, 
given infinite accuracy one could calculate the whole spectrum from one such correlation 
function. (Apart from states with accidentally zero overlap.) In the real world of finite 
accuracy, however, extracting more than the ground state energy, E , is an ill-conditioned 
problem so we shall proceed by the standard variational calculation described later on. 

To obtain Eq from a numerical calculation of C(t), we perform a fit with a single expo- 
nential to C(t) for t > t\ and choose for t\ the lowest value for which we obtain an acceptable 
fit. This is a way to estimate the asymptotic exponential decay in eqn([2]). Now, since the 
fluctuations that determine the error in the Monte Carlo calculation of C(t) are themselves 
proportional to a higher order correlation function which can be seen to possess a disconnected 
piece, the error is approximately independent of t, whereas C(t), as we see from eqn(j2J), de- 
creases at least exponentially with t. Thus the error/signal ratio increases exponentially with 
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t and one can only hope to obtain E accurately if the corresponding ground state already 
dominates C(t) at small values of t, i.e. one needs operators <ft that have a large projection onto 
the desired ground state. This can be achieved using standard smearing/blocking techniques 
(see e.g. [6J, [T] ) which produce operators that are smooth on physical length scales, with little 
computational effort. If clEq is small enough that the error/signal ratio is small for a nontrivial 
range of t = an t > t%, then the success of a single exponential fit provides significant evidence 
for it being appropriate. If however aE Q is so large that this range is small, then we lose such 
evidence, and typically we will be using a single exponential fit over a range of t where the 
contribution of higher excited states may still be significant. Thus our estimate of aEo will 
tend to become too high as aEo increases. This is an important sytematic error that will be 
visible in some of our results - in particular when we consider very long flux tubes (which are 
very massive) or highly excited flux tubes, even when they are not very long. 

To calculate the excited states as well as the ground states, we carry out the following 
variational calculation. Consider a set of n different smeared/blocked operators, {<fii(t) ; i = 
0,...,n — 1}, with the desired quantum numbers, and which have been normalised (but are 
not in general orthogonal). Now find the normalised linear combination of these operators, 
(pit) = $ (t), that maximises C{t = t ) = (t ) (f)(0)) = (0 f (O)| exp{-ift o }|</>(0)) where t 
is chosen small enough that all the correlators are accurately determined. (In practice we 
choose to = a.) Then <3>o provides our best estimate of the ground state wavefunctional. Now 
consider that part of the basis {4>i} that is orthogonal to <3>o ; and find the operator (j)(t) = $i(t) 
that once again maximises C(t = to) = (<^(0)| exp{— Hto}\(p(0)) but now over this reduced 
basis. Then $! provides our best estimate of the first excited wavefunctional. Repeating 
this procedure we obtain a set of operators, ; i = 0,1,...}, that are our variational 
approximation to the exact energy wavefunctionals, {^i}- For each <3>j we now consider the 
correlation function ($|(t)$j(0)), perform single exponential fits as described above for the 
ground state, and extract an estimate of aEi. 

As for the ground state, it is critical that one has good overlaps, and in practice that 
means that we need |(\l/j(0)$j(0))| 2 > 0.9. To achieve this for several excited states we need a 
large basis of operators. In practice we use 0(200) operators, as described below. For excited 
states there are additional systematic errors. First, the excited state may be unstable, e.g. if 
(Ei — Eo) is greater than the lightest scalar glueball mass. This will not lead to any visible 
effects in our calculation as long as the decay width is small. This will be the case at large 
N, where decay widths vanish, but even in SU(3) they are usually too small to be important. 
(Hence the difficulty in observing string breaking directly in lattice QCD.) Another systematic 
error arises from the fact that an excited trial operator, $j, may have a non-zero overlap onto 
a lighter eigenstate, i.e. |(\l/j(0) ( l ) j(0))| 2 7^ for j < i. If this overlap is small then it is only at 
large t that ($j(t)$j(0)) will behave like cx exp{— Ejt} and there will be an intermediate range 
of t where it will behave like oc exp{— Erf}, and where one can extract Ei. However as this 
overlap becomes more significant, the ambiguity in extracting Ej obviously will become larger. 
Another systematic error arises from the fact that the excited states show near-degeneracies 
that can be distorted by mixing effects which can be enhanced and complicated by the effects 
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of statistical fluctuations in our variational calculation. [J For all these reasons we shall confine 
ourselves to discussing only the lightest few excited states. 

2.3 winding flux loops 

Consider a loop of fundamental flux that winds once around the x-torus, so that it is of length 
I = aL x . A generic operator (pi that couples to such a periodic flux loop is the trace of an 
ordered product of link matrices, l p , along a space-like curve p that winds once around the x- 
torus, with the matrices in the fundamental representation. (This is an example of a Polyakov 
loop, or Wilson line, which may be in a general representation.) If we multiply all the U x link 
matrices in some given x-slice of the 3-volume by an element z G Zn, then l p — > zl p . This 
center symmetry keeps track of the net winding: a /c-string will transform as l p ^ — > z l Pt k- We 
will now specialise to the k = 2 strings which are the subject of this paper. 

We will consider k = 2 flux tubes of two kinds. First there is a flux tube that winds 
once around the torus but contains flux in a representation of SU(iV) that arises from the 
product of 2 fundamentals and any number of adjoints. The trace of a Wilson line in such a 
representation can in general be expressed in terms of a sum of operators that involve traces 
of powers and powers of traces of l p and of l p , such that the net power of l p minus the net 
power of l p is 2. In practice we shall limit ourselves to the operators Tr/ 2 and {Tr/ p } 2 whose 
linear combinations correspond to the totally antisymmetric (k = 2A) and totally symmetric 
(k = 2S) representations 

$ 2A = Trl 2 p - {Tr/ P } 2 ; $ 25 = Tr/ 2 + {Tr/ P } 2 (3) 

Since a k = 2S source can become a k = 2A source through gluon screening, there is no 
reason for such operators to be orthogonal or to assume that such SU(iV) representations will 
be useful in labelling the actual eigenstates. Whether they are useful is a dynamical question 
that we shall determine in this paper. The second type of k = 2 flux tube that we consider is 
one that carries fundamental flux and winds twice around the torus. A typical trial operator 
for such a flux tube will be 

$«,= 2 = Tr{l p l p ,} (4) 

where the two single-winding paths p and p' are in general not the same but the joint path 
p = pp' is continuous and periodic with a period 2l x . Note that the w = 2 operators need 
not be orthogonal to the k = 2A and k = 2S operators. Indeed because of the strong k = 2 
binding, one might expect that if p and p' are similar, so that the second winding of the flux 
tube is physically close to the first winding, the attractive interactions will give a significant 
probability of binding into a singly-wound k = 2 flux tube, and this will be reflected in the 
overlaps of the corresponding operators. One has therefore the heuristic expectation that to 
see a genuine w = 2 flux tube one needs to use very different paths p and p' . 

In order to obtain good overlaps onto the true energy eigenstates, we use smeared/blocked 
links in our construction of l p . Usually, and always when it matters, as in the projection onto 

1 For example, a correlation between two nearly-degenerate states, that would average to zero with infinite 
statistics, can induce an artificial level repulsion. 
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specific representations of the flux, we make sure that our blocked matrices are all SU(iV). 
(For historic reasons this is not always the case for some of our oldest calculations, although 
this will not affect any conclusions we shall draw from them.) Using operators with different 
levels of smearing probes different transverse size scales, and this forms part of our variational 
basis. 

What are the other quantum numbers that label such a flux tube? First, we recall that in 2 
space dimensions, there are no rotations around the axis of the tube and in the confining phase 
a flux tube that is rotated by n/2, so that it winds around the y-torus, will have no overlap 
with one winding around the x-torus, so this symmetry is uninteresting. We do however have 
translations and corresponding momenta. We project onto zero-momentum transverse to the 
flux tube, p y = 0, by adding all y-translations of our basic operator. We do not consider p y ^ 
since previous calculations have shown that for the range of lattice spacings we shall consider, 
the continuum energy-momentum dispersion relation is accurately satisfied and so there does 
not seem to be anything new that one might learn using non-zero transverse momenta. More 
interesting is the momentum along the flux tube, p x = 2ixq/l. Here we translate the loop by Xq 
in the x direction, multiply it by exp{ip x xo}, and sum over all such translations. Clearly if the 
flux tube is invariant under translations along its axis, as we might expect if it is in its ground 
state, then imposing q ^ in this way will produce a null state. To have q ^ the state must 
be excited in some way. Thus calculating the spectrum for various q ^ promises to teach us 
something new. Another symmetry is (two dimensional) parity, P. Again one expects that to 
have P = — one requires a deformation of the flux tube, i.e. an 'excitation', that is different 
under reflection. Finally there is charge conjugation C . This operation changes the direction 
of the arrow, l p — > l^ p , so that Re0 is C = + and Im0 is C = — . Since l p — > zl p and Z| — > zU^ p 
under the center transformation introduced earlier, there is in general zero overlap between l p 
and l p so the C = + and C = — sectors are degenerate and the label is therefore uninteresting. 
However this is not the case when k = N/2. In that case there can be mixing and the C — ± 
sectors need not be degenerate. In our case, that occurs for k = 2 in SU(4). However it is 
easy to see that, for SU(A), Im$ fc=2 A = 0. Thus to the extent that there is a k = 2A sector 
of states, we will only have C = + therein. For k = 2S on the other hand, C — ± can be a 
useful label in SU(4). 

In Appendix A we describe in more detail the construction of the actual operators we use. 

3 Flux tubes as strings 

The immediate goal of our calculations is to provide information about the spectrum of k- 
strings that will complement what we have learned previously about fundamental k — 1 flux 
tubes. Our wider purpose, however, is to provide useful information for theoretical approaches 
to understanding confinement and the dynamics of QCD. Thus a central step in our analysis 
is to compare our results to what one might expect on very general grounds, so as to pin-point 
what might be new and interesting. In this section we shall first discuss what these general 
expectations might be, and then we summarise the detailed spectrum of the free Nambu-Goto 
string theory (in flat space-time) since this turns out to provide a very useful comparative 
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benchmark. 

The flux tubes that will interest us here, are those that are wound around a spatial x-torus 
with all other space-time dimensions effectively infinite. As we shrink the length of this torus, 
I = aL x , such a system will undergo a phase transition at a critical value I — l c — 1/T C where 
T c is the deconfining temperature. For I < l c the ground state is no longer a flux tube, but 
rather has zero energy. In D = 2 + 1 and for N > 4 one finds l c y/&f — 1-1 |8J. Thus if we 
want to discuss closed flux tubes, we must choose I > l c . 



3.1 general expectations 

Consider a flux tube of length I = aL x winding around the x-torus. Translations leave its 
energy unchanged. However under transverse translations the state is not the same: the 
choice of a position for the flux tube breaks this symmetry spontaneously and we expect to 
have corresponding massless fluctuations in the excitation spectrum of the flux tube. Since 
we would expect the ground state to be invariant under translations along the flux tube, we 
do not expect this longitudinal translation symmetry to be spontaneously broken and we do 
not expect any further Goldstone bosons from that source. 

Now, since the SU(iV) field theory has only one scale, say l a = 1/ y/oj, the width of the flux 
tube should in general be on that scale. On length scales ^> l a , and for excitation energies 
<C l/l a the flux tube should effectively behave as a thin string. The classical fluctuations 
of such a periodic string have wavelength A n = l/n and hence frequency u n = 2im/l. The 
corresponding oscillators will be quantised, and the excitations can be thought of as made 
up out of phonons of energy uo n and with momenta p n = ±uj n along the string. This is 
the expected Goldstone boson, which is massless and can have any momentum allowed by 
periodicity. At large / all these eigenstates converge 

lim E n /al = 1, Vn. (5) 

At finite / there will be shift in the ground-state energy from the zero-point energies of the 
oscillators (essentially a Casimir energy) and this leads to the well-known Luscher correction 
[9] to the string energy: 

which reflects the asymptotic low-lying n = excitation spectrum 

The coefficient of the 0(1/ al 2 ) term is universal and the value shown in eqn (!6|7|) is the one 
for the simple bosonic string universality class, where the only massless excitations of the flux 
tube arise from the spontaneously broken transverse translation invariance. There is strong 
numerical support for this choice in the case of both fundamental and k = 2 flux tubes [HE], 
and this is of interest since in principle there might have been other, less obvious, massless 
fluctuations that would have led to a different coefficient. 
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The large-Z expansion in powers of 1/er/ 2 , that we see in eqns fl6|7|) . arises naturally in the 
classic framework of Luscher, Symanzik and Weisz, [9j[T0] where one categorises and constrains 
the possible interaction terms in the transverse displacement fields. These interaction terms 
will in general lead to non-trivial interactions amongst the phonons and are encoded in the 
higher order corrections to the ground and excited state energies shown in eqns( l6|71) . This ap- 
proach has recently led to the remarkable conclusion [10] that the next term in this expansion 
is oc 1/(<t/ 2 ) 2 with the same coefficient as in the Nambu-Goto spectrum (see below). A simi- 
lar conclusion was independently reached in [11] using the alternative Polchinski-Strominger 
framework [12] where one builds up an effective D-dimensional string theory within a general 
string path-integral setting. 

These two theoretical approaches have in common that they only hope to treat very long 
flux tubes and, in effect, only eigenstates that are string-like in the sense of eqn([5]). However 
it is common to think of the confining flux tube as being qualitatively similar to a Nielsen- 
Olesen vortex, albeit one that is a non-Abelian dual-superconducter. Such a vortex will 
in general also have a pattern of excitations that reflects the effective theory producing its 
non-trivial structure. (For an explicit example in related, more tractable theories see [13J.) 
These excitations will typically be massive, leading to string eigenstates with energies that 
will typically be 'much' larger than those of stringy excitations at large enough /: 

l^oo An f D — 2\ 

E(l)~E (l)+n > E n (l)~al + T ln-^ r )+... (8) 

One's first guess for /i might be something on the order of the lightest glueball mass, or the 
splitting between that state and its first excitation. Using values for the lightest glueball 
masses from [6] this would suggest 

H ~ 2 - 4y/a f . (9) 

Such states should become particularly visible at smaller values of I where the gap between 
the stringy states becomes larger than this estimate for \x. 

One might simply give up on describing such massive excitations in the belief that the spec- 
trum of anything but a very long flux tube will be too complex to be theoretically tractable. 
However our earlier calculations of the low-lying excitations of the fundamental flux loop have 
revealed that the simple stringy predictions of the free bosonic Nambu-Goto string theory 
hold very accurately even for flux tubes as short as ly/aj — 1.5 [TJ [2]. That a flux tube which 
is only a little longer than it is wide should display simple string-like excitation modes is not 
something that one can understand within a generic vortex picture. It suggests that despite 
appearances short flux tubes do have a stringy description. Does this mean that the flux tube 
is indeed a thin string, with its apparent width simply a manifestation of the zero-point fluctu- 
ations described above? This would immediately come up against the difficulties of defining a 
consistent bosonic string theory in 3 or 4 dimensions. A more interesting possibility is offered 
by a dual gravity description of the AdS/CFT type. (In fact non-CFT and possibly non-AdS.) 
Here the flux loop is indeed always described by a string but now this string is in a geometry 
that is, at least initially, 10 dimensional. Of course the spectrum is now that of a Nambu-Goto 
string in curved rather than flat space-time, and this curvature is crucial to providing linear 



8 



confinement rather than a conformal-like Coulomb interaction. In this picture it may not be 
unnatural for some stringy states to arise even for short strings, and the massive modes are 
of great interest as they provide clues as to the geometry of the dual theory. 

3.2 Nambu-Goto closed string spectrum 

The simplest specific bosonic string theory is the free string theory: Nambu-Goto in flat space- 
time. Of course such a theory can only be consistently formulated in 26 dimensions. However, 
as discussed by Polchinski and Strominger [12] , if one is interested in an effective string theory 
for single long strings, then the Weyl anomaly does not pose an obstacle to working in D < 26. 
(As was already pointed out long ago by Olesen [Hj.) The spectrum of the Nambu-Goto free 
bosonic string theory inD = 2 + lorD = 3 + l and in the sector of a single long string, has 
also long been known [151 HE]- We shall briefly summarise this spectrum here, since it will 
provide our central point of comparison in this paper. 

Consider a single closed string that wraps w times around a spatial torus. The free 
string spectrum corresponds to quantised oscillations of the string that may be thought of 
as phonon modes traveling clockwise (right-handed) and anticlockwise (left-handed) along 
the background string. Our convention is to take the momenta of the latter as positive and 
of the former as negative. The modes are massless, corresponding to the Goldstone bosons 
that arise from the fact that the presence of the background string spontaneously breaks the 
transverse translation invariance. So the energy of a phonon equals the absolute value of its 
momentum. This momentum is along the string and is constrained by the latter's periodicity 
to be ±2nk/lw where k = 0, 1, 2, . . . and Iw = aL x w is the length of the string. Of course, 
the total momentum is constrained to be 2irq/l = 2ixqw/lw, with q = 0, ±1, ±2, . . . , since the 
Hamiltonian is defined on a spatial torus of size / (and not Iw). 

Let nL{k),nn(k') be the number of left and right handed phonons with momentum 27rk/lw 
and —2ixk'/lw respectively. A free string state is defined by the total set of phonons and the 
total energy and momentum will be determined by 

k n L (k) k' n R (k') 

The sum of the phonon momenta is just the total momentum of the excited string and is 
therefore constrained to be 

N L -N R = qw. (11) 

In addition the states have a parity, P. Under P (in two space dimensions) the transverse 
displacement changes sign, so the corresponding creation operators for these displacements, 
the right-handed and left-handed phonon operators, will change sign. That is to say, 

p = /-Qnmnber of phonons _ (_]\12k Y.n L (k) n L (k)+J2k' En fl (t') ^2) 

Writing down the 4-momentum of such a state and taking its scalar product, one obtains the 
energy-squared: 

p2 / / x2, B (Nl + Nr D-2 \ , {27rq\ 2 

E N L ,N R , q , w = {<rlw) + 8na\^ —j + J (13) 
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where the negative piece in the middle term arises from the zero-point energies of the oscilla- 
tors. These energy levels will have a degeneracy that is determined by the number of different 
ways one can produce the same value of Nl + Nr in eqn( fT0l) . As an example we present in 
Table [1] the seven lowest energy levels of a closed bosonic string that winds once around the 
torus for q = 0,1, 2, in terms of the left and right movers (creation operators) and a—y. 

3.3 corrections to the free string 

We now consider how to incorporate corrections to the free string spectrum. We consider 
q = and w = 1 but the extension to other values should be obvious. 

The conventional Luscher analysis [TO] of confining strings focuses on the effective La- 
grangian for the displacement field, and categorises the interactions in powers of these fields 
and derivatives. Originally this showed that in an expansion of E{f) in powers of 1/al 2 , the 
first correction, which arises from the oscillator zero-point energies, is universal and is propor- 
tional to the string central charge. Given the universality it is no surprise that this term, the 
well-known Luscher correction, is precisely what one obtains as a first correction to al when 
one expands the Nambu-Goto expression for E(l) in eqn f)13l) in powers of 1/al 2 . More recently 
[TO] this analysis has been carried further, and it has been shown, on quite general grounds, 
that in D = 2 + 1 the next term in this expansion also coincides with what one obtains from 
Nambu-Goto, i.e. 



where n = 0,1,2.... Quite independently, an extension [llj of the work of Polchinski and 
Strominger [T2] reached the same conclusion (but for D = 3 + 1 as well) . 

In our calculations of the k = 1 flux loop [TJ [2] we discovered that the actual spectrum 
is reproduced by the free string Nambu-Goto spectrum much better than one would expect 
from eqn (fl~4"]) . In fact there is even good agreement for values of / that are so small that the 
expansion of eqn (fl~3]) in powers of 1/al 2 completely diverges. The implication is clear [2]: we 
should use the full Nambu-Goto expression as our starting point, and consider corrections to 
that. In that spirit we write 



where p = 3 corresponds to the 0(l/l 4 ) correction in eqn f[T4"]) . In this paper we will fit our 
k = 2 spectra with eqn (!15p . with the leading permitted correction of p = 3. 

4 Spectrum of k = 2 flux tubes 

We have performed calculations in SU(4) and SU(5) at j3 = 50 and (3 = 80 respectively. Since 
(3 = 2N 2 /X, where A = g 2 N is the 't Hooft coupling, these values of (3 should correspond to 




(14) 




P > 3 



(15) 
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roughly the same value of the lattice spacing a, and indeed we find 




0.1310(2) 
0.1298(1) 



SU{4), (3 = 50 
SU{5), p = 80 



(16) 



Extensive experience with glueball masses [6] the deconfinement temperature [8j and the 
spectrum of fundamental strings [U [2] indicates that lattice spacing corrections for such values 
of a are negligible and we shall therefore assume that what we obtain is indeed continuum 
physics. 

We calculate the k = 2 spectra of flux tubes winding around a;-tori that range in length 
from I = 12a to I = 32a. As we see from eqn fU6l) this corresponds to l\f<Jf € [1.56, 4.20]. Thus 
our shortest flux tubes are presumably not much longer than they are wide - hardly 'tubes' 
at all. For comparison, the minimum length l c below which we lose such flux tubes completely 



Since we are interested in accurately identifying corrections to the asymptotic linear variation 
with I, it is important that we minimise systematic errors that vary with I. Thus at small 
I we must make L t large enough so that the contribution of both k = 1 and k = 2 loops to 
the partition function is negligible, and all our energies are normalised to the correct vacuum 
energy. Similarly we must also make the transverse size L y large enough to avoid finite size 
corrections. Thus, following a variety of checks, the lattices we use for SU(5) are: 12 x 32 x 80, 
16 x 32 x 40, 20 x 24 x 32, 24 x 24 x 32, 28 x 28 x 32, 32 x 32 x 32. For our SU(4) calculations, 
which were performed earlier, we use slightly smaller lattices, where finite volume corrections, 
for the smallest values of / could be non-negligible: 12 x 24 x 40, 16 x 24 x 40, 20 x 24 x 32, 
24 x 24 x 32, 28 x 28 x 32, 32 x 32 x 32. We typically perform 500,000 Monte Carlo sweeps 
and measure the full correlation matrix every 5 sweeps. 

As described earlier, the relevant quantum numbers in labelling our k = 2 flux tubes are 
the momentum, q, along the flux tube and the parity P. We ignore charge conjugation since 
it is only relevant for N = A, and then only for the heavier k = 2S states. In addition, the 
flux may be in different representations of SU(iV). The SU(N) representation is, however, 
not a conserved quantum number since it can change under screening by gluons. Nonetheless 
we shall see later that the mixing between k = 2A and k = 2S states is typically very small, 
and so in our analysis in this section we shall categorise the states in this way as well. (Our 
basis of operators does not allow us to explore other representations.) Note that where it 
significantly improves the calculation we use our full k = 2 basis of operators, although we 
categorise the states by their dominant k = 2A or k = 2S component. 



In Figs. [T] and [2] we plot the energies of the lightest states with momenta q = 0, 1, 2 along the 
flux tube, for SU(4) and SU(5) respectively. We observe that for all values of I the lightest 
q = state has P = +, the lightest q = 1 state has P = — , and for q = 2 there are two 



is [8] 




8.1a 
8.3a 



SU(A), (3 = 50 
SU(5), [3 = 80 



(17) 



4.1 k = 2A 
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nearly degenerate lightest states, one with P = + and one with P = — . We note that these 
quantum numbers and (near-) degeneracies are precisely those of the Nambu-Goto free string 
theory. (See Table [U) All these states have a projection onto 2 A that is very much larger 
than onto 2S and we therefore regard them as constituting part of the 2A spectrum. We have 
removed all lattice units by normalising the loop energies and the loop lengths to the value of 
(i\/crf that has been obtained from k = 1 flux loops in the same calculation. (The errors on 
this quantity are much smaller than all our other errors.) 

We fit the q = ground state using the Nambu-Goto expression with the leading permitted 
correction, as in eqn (fT5l) . [TU1 fTT] . 

E ^=^i i -i^)-^j^=v (i8 » 

to obtain 

a 2 a 2A = 0.023233(51) = 1.354(5)aV / ; C 3 = 2.93(21) : SU(A) (19) 

a 2 a 2 A = 0.025825(50) = 1.533(4)aV / ; C 3 = 2.28(22) : 517(5) (20) 

Note that the values of o^a/* 7 / are close to the quadratic Casimir ratios, k(N — k)/(N — 1). 

We observe that when expressed in terms of natural units, as in eqn (fT8|) . the correction 
to Nambu-Goto has a coefficient C 3 ~ 0(1). This is a natural value if we regard the string 
theory as an effective theory, in the Nambu-Goto universality class. It is to be contrasted with 
the value C% ~ 0(1/10) that one finds for the fundamental flux tube, and serves to emphasise 
how remarkably well the latter is described by the free bosonic Nambu-Goto string theory. 

Using these calculated values of <T2A, we calculate the lightest q = 0,1,2 energies of the 
Nambu-Goto model, and plot the results in Figs [1] and [2j We observe an excellent agreement 
between our numerical results and the Nambu-Goto values. Note that once the q = ground 
state has been fitted, the q = 1, 2 energies are parameter- free predictions. 

Since the Nambu-Goto prediction appears to work so well, it is worth examining the results 
more closely to try and isolate how well the phonon excitation energy is being reproduced by 
our numerical results. To do this we note from eqn (fT3|) that the piece of the energy that is 
due to the phonons is given by 

A(q, I) = E 2 gs (q; I) - E%{q = 0; I) - {^f (21) 

Using our numerically calculated values for E gs (q; I) we obtain the results displayed in Figs. [3] 
and H]for SU(4) and SU(5) respectively. (We choose to use our numerical value for E gs (q = 0; I) 
rather than the Nambu-Goto prediction since the former already has some small corrections 
which we would regard as belonging to the Casimir energy rather than to the excitation en- 
ergy that we are trying to isolate here.) For comparison we show the Nambu-Goto phonon 
excitation energies as well. We observe that the simple free string excitation energies pro- 
vide an excellent description of the string excitation energies over almost our whole range 
of string lengths, with significant deviations only beginning to appear at the smallest value 
of I, ly/a ~ 1.5. Given the qualitative (quantum numbers, degeneracies) and quantitative 
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agreement between Nambu-Goto and our results, we will assume that these eigenstates are 
very close to those of the free string theory, with the corresponding phonon excitations. 

We began this Section with the ground states for various q because, as explained earlier, 
ground state energies can be extracted more reliably than excited state energies. We now turn 
to the excited states that form the low- lying q — 0, P = + spectrum. We plot our results for 
SU(5) in Fig. [6] together with the Nambu-Goto predictions. (The SU(4) results are essentially 
the same and are given in Fig. [5j) The first excited state shows substantial deviations from 
the free string energies, and the next two states even more so. However at least the former 
does seem to be approaching the free string energy as I increases. This might tempt us to 
see it as a string excitation, but this would be dangerous, since we also observe that to quite 
a good approximation this first excited state is greater than the ground state by a constant 
value, 

E 1 (q = 0,l)~E o (q = 0,l)+e ; e~2^a f , (22) 

which is reminiscent of our expectation for non-stringy states as described in eqns (181191 . Such 
a state would cross the level of the first string excitation at some larger I, and perhaps that, 
not a convergence, is what we are seeing in Fig. [6j At this stage we clearly cannot decide how 
to classify this excited state. We shall return to this question in a later Section, where we 
shall find that looking directly at the wavefunctional will give us an unambiguous answer. 

4.2 k = 2S 

As we shall see later on, all the states we analyse appear to transform primarily as k = 2A 
or as k = 2S (or as w = 2 fundamental). However the lowest k = 2 A states are much lighter 
than the corresponding k = 2S states, so the latter will be much more susceptible to the 
systematic errors discussed in Section I2T21 - in particular to their masses being over-estimated. 

Bearing this caveat in mind, we display in Figs. Hand [H]the energies of the lightest k = 2S 
states with momenta q = 0,1, 2, for SU(4) and SU(5) respectively. Qualitatively the picture 
is much as for the k = 2A states: we find that for all values of I the lightest q = state has 
P = +, the lightest q = 1 state has P = — , and for q = 2 there are two nearly degenerate 
lightest states, one with P = + and one with P = — . All this fits in with the Nambu-Goto 
free string picture. Turning to a quantitative comparison with Nambu-Goto, we encounter an 
immediate difficulty. Our estimate of the / = 32 ground state energy is too high to fit in with 
an E(l) tending to an asymptotically linear behaviour. We therefore assume that what we are 
seeing here is the result of a systematic error of the kind described in Section 12. 2\ and which 
is perhaps twice the size of the statistical error. If we exclude this value from our fit to the 
q = ground state using eqn (TT3"j) . we obtain 

a 2 a 2 s = 0.03957(45) = 2.306(27)aV / : S77(4) (23) 
a 2 a 2S = 0.03793(33) = 2.251(20)aV / : SU(5) (24) 

(where the correction coefficients, C3, are not well-determined). In fact if we had included the 
I = 32 data point we would have obtained almost identical values of a 2 s, but with a much 
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worse x 2 - We note once again that the values of cr2s/ cr f are close to the quadratic Casimir 
ratios, in this case k(N + k)/ (N + 1). 

Using these values of a^s, we show the corresponding q = 0,1,2 ground state Nambu-Goto 
spectra in Figs. [7] and El Semi-quantitatively the agreement is good; however the longer, 
more massive states tend to overshoot the predicted values. It is plausible that this is again 
due to the systematic errors of our calculation growing with the calculated mass. 

In Fig. [9] we extract the phonon excitation energies from our SU(5) values, using eqn (!22|) . 
(The SU(4) results are very similar.) There is good qualitative agreement with Nambu-Goto, 
but clearly we cannot say more than that. 

In contrast to the k = 2A case, we have severe difficulty in extracting excited state energies 
and hence do not attempt to produce a plot that is analogous to Fig. [6j 

To conclude, it is probably fair to say that the q — 0,1,2 ground states provide significant 
evidence that the k = 2S flux tube also behaves like a free bosonic string, at least for these 
states, but that any attempt at a stronger statement is hampered by the current uncertainties 
of our numerical calculations for these very massive excited states. 



5 Group or centre? 

For each k, and for each set of conserved quantum numbers P, C and q, there should only 
be one absolutely stable flux tube. (Except e.g. where it can decay into a string with the 
same k but with some of the other quantum numbers different, and accompanied by glueballs 
that carry the compensating quantum numbers.) However, flux tubes carrying flux in an 
SU(N) representation with iV-ality k that is, in general, not stable, can be stable enough to 
be readily identifiable in our calculations. If so it is an interesting question to ask if they are 
individually described by an effective bosonic string theory with an appropriate string tension. 
We have, of course, already answered this question in the affirmative in earlier sections, where 
we examined the k = 2 flux tubes in the totally symmetric and antisymmetric representations. 
(Which are the only ones accessible to our limited choice of operators.) Here we will present 
evidence that the portion of the k = 2 spectrum accessible to us is, to a good approximation, 
the sum of disjoint k = 2A and k = 2S spectra. This provides the promised continuation of 
the discussion of this question in [5] to which we refer the reader for further background. 



5.1 operator overlaps 

We begin by looking at the overlaps of our k = 2A basis operators, {$ 2j 4> * = 0, 1, 2, ...}, onto 
their k = 2S counterparts, {$2,s^ = 0,1,2,...}. For this purpose we define the normalised 
overlaps: 

n r ■ ,x < $Lu(*) WQ) > , oK . 

OAS(l,J',t) = 1 ! — ; 1 (25) 

< $i Ai (m 2A>i (o) >V2< ^ 2S)j ( t )^ SJ (o) >V2 

In Table 4 of [5] we presented results for such overlaps, at t = 0, using a small basis of 
blocked/smeared Polyakov loops of length / = 24a on a 24 2 32 lattice in SU(5) at (3 = 80 
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(corresponding to ly/vj ~ 3). We found that most of the overlaps were zero within errors, 
except for those where at least one of the operators was at the largest blocking level; but even 
there the largest overlap was still a very small ~ 0.03. This provided some striking evidence 
that the theory does indeed break up into almost disjoint k = 2A and k = 2S sectors. However, 
while this small basis has a reasonably good overlap onto the k = 2A and k = 2S ground 
states, the overlap onto excited states is in general poor. It is therefore interesting to see what 
happens when we consider the much larger operator basis of the present paper, which has a 
good overlap onto most of the low-lying states. Another question is how much the results of 
|5| are affected by finite volume corrections. The point being that the only significant overlaps 
were for the largest 5'th blocking level, where the blocked link is of length 16a and has a width 
that is at least as large - that is to say, the operator more-or-less covers the whole spatial 
torus. Finally a very interesting question is what happens for t ^ 0, since this probes more 
directly the represntation content of the energy eigenstates. 

To address the question of finite volume corrections we perform some t — calculations in 
SU(4) at (3 = 32. (Since in SU(4) the k = 2A operators are necessarily C = +, we take their 
overlaps with the C = + piece of the k = 2S operators.) We start with / = 16a because this 
also corresponds to ly/oj ~ 3. Note that on a 16 2 spatial lattice the 5'th blocking level really 
does cover the whole spatial torus; so any finite size corrections should be larger than in [5]. 
We show a selection of the overlaps on a 16 2 20 lattice in Table [21 including all those which 
are large enough to be significant. (Nearly all those not shown are zero within errors.) We see 
that a few overlaps are indeed not small, but that they again occur when one of the operators 
is at the largest blocking level. We now see what happens if we keep the same loop length, 
but vary the transverse spatial size, i.e. we work on 16 x l y spatial lattices where we vary l y . 
What one finds is that the overlaps rapidly decrease with increasing l y (and rapidly increases 
as we decrease l y below Z).We show examples in Table [2] for l y = 20 and l y = 48. The decrease 
is however not to zero but to a finite value; in fact the l y = 48 values are asymptotic within 
errors. These values are very small, and much smaller than on the original 16 2 torus. We now 
turn to what happens as we increase I while staying on an I 2 lattice. As we see in Table [2] 
this also strongly decreases the overlaps. In this case all the overlaps appear to tend to zero 
as I increases. Of course higher blocking levels will become possible at larger /, and we would 
expect the largest of these to sometimes generate significant overlaps, but such operators will, 
eventually, have very small overlaps onto the low-lying flux loop states. 

We now turn to the overlap between all the operators in our extended basis. These have 
been calculated on a 32 3 lattice for P = +, q = and SU(4) at (3 — 50. A Table is not 
practical, so we provide in Fig. [TO] a three-dimensional plot of the values of \OAs(hj',t — 0)1 
as i and j run over our basis of some 80 operators. As we might have expected from our 
more limited study above, the overlaps are mostly very small, and usually zero within the 
errors. Only those that involve operators at the very largest blocking levels are significant 
although even these are numerically small. (Here the 5'th blocking level is the largest, and it 
appears in our operator labeling at multiples of 5.) We can further assume, given what we 
found above, that most of these larger values are much enhanced by finite volume corrections. 
Thus this study confirms that even when we extend the basis so that it is quasi-complete 
for the low-lying spectrum of the k = 2 flux tube, the near-orthogonality of the k = 2A and 
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k = 2S sectors is maintained. Physically this suggests that the effects of screening k — 2S 
down to k — 2 A are very weak, and why that should be so, even in SU(4), needs to be better 
understood. 

Of course the near orthogonality of the k = 2A and k = 2S sectors does not necessarily im- 
ply that the energy eigenstates belong to one sector or the other (although it is a condition for 
that to be possible). To be sensitive to the eigenstate content we need to look at \C>As(hj] t)\ 
for t ^ 0. Since the correlators decrease exponentially with t, we only have accurate results 
for the lowest values of t and we will therefore restrict ourselves in the following dicussion to 
t < 2. Now, if the eigenstates were linear combinations of comparable amounts of k — 2 A 
and k = 2S, then we would expect to find \OAs(hj',t 7^ 0)| ~ 0(1). If on the other hand 
\OAs{h j',t 7^ 0)| is 'small' then this suggests that the lightest states are largely in one sector 
or the other. To determine which is the case, we turn to our high statistics calculation on 
a 24 2 32 lattice in SU(5) at (3 = 80. In Table [3] we show the same selection of overlaps as in 
Tabled but now for t — 0, 1, 2. (The t = results differ from those shown in [5] because the 
statistics here are ten times greater.) We see that the t ^ overlaps do indeed remain very 
small. Moreover some lower statistics calculations on a 24 x 32 x 32 lattice show that with such 
a larger transverse size, the overlaps decrease by a factor ~ 5 — 10. Thus we conclude that 
the overlaps for t ^ are also very small, and the low-lying eigenstates cannot be comparable 
mixtures of k = 2A and k = 2S. To make this more precise we shall now turn to an analysis 
of the actual eigenstates. 

5.2 state overlaps 

What we will do in this Section is to estimate the overlaps of the actual string eigenstates onto 
the k = 2A and k = 2S sectors. Our approach to this question, which is somewhat different 
to that employed in [5], is as follows. 

Let {^i} be the true energy eigenfunctionals ordered, as usual, by their energy {Ei}. 
Using only our k = 2 A basis, let the variationally calculated trial eigenfunctionals be {$2A,;}- 
Suppose now that <&2Aj corresponds to where j < i. (The possibility j < i arises when 
there is an eigenstate \l/ with a small enough overlap onto our k = 2A basis, that it does not 
lead to a corresponding $2A, so that the set of variational trial eigenfunctions is incomplete.) 
Then the normalised correlation function of $2Aj will be ~ exp{— E{t} in some range of 
t - the effective energy plateau. Then 7 2 a gives us the overlap of the eigenstate on $2A,j 
and hence, to a good approximation, on our k = 2A basis. Doing the same with the k = 2S 
basis should give us the corresponding overlap 72s. Now returning to reality, we know that 
we will not see such a state at all in any sector where its overlap is small. So in practice we 
can usually only calculate either 72A or 72s this way, but not both. So suppose it is 72A that 
is large enough to be clearly calculable. Then instead of looking within the k = 2S basis to 
estimate 72s, we look at the joint basis {k = 2A} ®{£; = 2S}. Here the overlap is 72AS and 
must be larger than 72A, so the state will certainly be visible to the variational calculation. If 
we take the k = 2A and k = 2S bases to be approximately orthogonal to each other, which we 
have shown above to be the case, then we can estimate — I2AS ~ i\a- ^ ^ ne state is more 
k = 2S than k = 2A then we follow the same procedure but with 2A with 2S interchanged. So 
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in this way we obtain an estimate of and 7I5 for each energy eigenstate. These estimates 
are of course approximate. In particular, our k = 2A and k = 2S bases are not complete, so 
that 72,4 + J2A < 1 an d becomes more so (typically) for higher excited states. In addition the 
non-orthogonality of the k = 2S and k = 2A sectors introduces an ambiguity in interpreting 
the results, but this is very small. 

We show a selection of such analyses in Fig. [TTJ The examples are all from SU(4) at 
j3 = 50. (The SU(5) results would be very similar.) We show results for the P = +,q = 
spectrum for the two string lengths shown, as well as results for the longer string length for 
the P = — ,q = and P — — ,q — 1 sectors. For each state we show our estimate of the 
k = 2A (red circle) and k = 2S (blue square) overlap squared. We see that all the states are 
predominantly in one sector or the other. In some cases the sum of the overlaps is significantly 
less than unity, and that simply reflects the overlap of that state on our whole basis. This 
effect becomes stronger for the higher excited states. 

6 cu = 2 'unbound' fundamental flux loops 

So far we have focused on k = 2 bound states. However the k = 2 sector should also include 
states that consist of a fundamental flux loop that winds twice around the spatial torus, and 
whose spectrum will be given, to some first approximation, by eqnf ll3l) with uj = 2. These 
states can be thought of as corresponding to two unbound fundamental flux-tubes. Such 
states might suffer large corrections because in D = 2 + 1 such a double winding loop will 
necessarily cross itself and we know that there are interactions between loops that are strong 
enough to bind them into the k = 2A strings. Moreover the u = 2 loops should typically be 
more massive than the corresponding k = 2A states and so they will typically appear in a 
part of the k = 2 spectrum where they might be confused with other excited states. All this 
should make them non-trivial to observe and so what we shall attempt in this Section is no 
more than a preliminary search for such uj = 2 flux loops. 

For this purpose we include in our basis additional operators of the form in eqn(JH). We then 
perform the same calculation as in the Section[4]but within this extended basis. Comparing the 
resulting spectra to those previously obtained in Section UJ we attempt to identify new states. 
The new basis of operators that has been used in our calculation consists of approximately 
~ 240 operators. 

In FigsJT2l and [131 we compare the S77(4), P = — , and q = spectra obtained with and 
without the extra uj = 2 operators, and we do indeed observe that the former appears to 
contain two extra states, which we have shown as solid black symbols. For comparison, the 
blue line represents the NG prediction for a uj = 2 k = 1 string, using a = <jf and Nr = Nl = 2 
in eqn fU3l) . The red line represents the NG prediction for the u = 1 k = 2 antisymmetric 
representation with the same values of Nl and Nr and with a = <Jia- As we see from Table UJ 
the ground state of a winding string with P = — and q = should be two-fold degenerate, 
and previous calculations for uj = 1 k = 1 flux tubes [2] have shown that what one gets is 
one state which is close to the theoretical NG prediction and another one that approaches 
the theoretical prediction from above as I increases. We observe precisely such a behavior in 
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Figure [121 with the lighter state, represented by a star, being close to the w = 2 prediction 
and the heavier state, represented by a circle, approaching the blue line from above. 

To confirm that these two states are primarily associated with the extra uo = 2 operators 
in our extended basis, we show in Figure [TH the projection in of the variational eigenstates 
corrsponding to the ground and first excited states onto all the operators in our basis. We 
have chosen to look at / = 16a (1-JoJ ~ 2.1) where these states are (supposedly) the lightest 
u = 1 k — 2A and u = 2 k = 1 states. For this lattice size we have 160 operators in total, 
of which the first 100 are our ordinary k = 2 operators, and the last 60 are our extra u = 2 
operators. As can be seen, the ground state projects mainly onto ordinary k = 2 operators, 
while the first excited state projects almost entirely onto the w = 2 sector of the operator 
basis. This supports our earlier interpretation of these states. 

It is interesting to note that the w = 2 operator that contributes the most to the first 
excited state is the rectangular pulse (Table HI number 3) at blocking level 4 and, therefore, the 
transverse deformation has a physical length of 16 lattice spacings. A pictorial representation 
of this operator is included on the panel of Figure dH We see that it is shaped like an oscillating 
closed string with two nodes, doubly wound around the torus. This provides further heuristic 
support for our conclusion that the complete k = 2 string spectrum includes unbound wave- 
like uj = 2 states. 

7 Looking for extra massive states 

As we remarked in Section 14. 1[ the interpretation of the first excited state in the k = 2A, 
P = +, q = spectrum shown in Fig. [6j is ambiguous. It might be the first string excitation, 
with a large correction to the Nambu-Goto free string energy, which is shown in Fig. [6] and 
which it appears to be approaching as I f. Alternatively it might be a new excitation that 
arises from one of the massive modes of the ground state flux tube being excited. Indeed, as 
we see from Fig. [6], the energy of this state does differ from that of the ground state by roughly 
a constant as I varies. In this scenario, the energy approaches the energy of the first excited 
Nambu-Goto state, only to cross it at some larger value of I and in that case it is presumably 
the next state up in energy which is the true stringy excitation (and which indeed is much 
closer to the Nambu-Goto value for most of our values of I). 

Since the identification of such extra states would (begin to) provide a completely new 
source of useful information about the string description of the gauge theory, resolving this 
ambiguity is important, and that is what we turn to now. 

If the state is indeed an approximate Nambu-Goto-like string excitation of the k = 2A 
flux tube, then we would expect its wave-functional to have the appropriate 'shape'. What 
that 'shape' should be, in terms of our highly blocked/smeared link matrices, is not at all 
evident, but it is something we do not need to know. It is enough to note that we have 
already established [2] that the relevant part of the fundamental k — 1 flux tube spectrum is 
Nambu-Goto-like. Indeed the k — 1 analogue of Fig. [H] shows only small corrections to the 
free string energies. So all we need to do is to compare the wavefunctionals of the first excited 
k = 1 and k = 2A states, in exactly the same lattice calculation, and see if they are very 
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similar or very different. 

The way we make this comparison is as follows. Let {(fif, i — 1, n } be our set of winding 
operators where the number is typically n Q ~ 100. In the present case, we choose them to have 
P = + and q = 0. These operators are group elements (not yet traced) and may be in any 
representation of SU(iV). Suppose the flux is in the representation 1Z. When we perform our 
variational calculation over this basis, we obtain a set of wavefunctionals, which are an 
approximation to the corresponding eigenfunctions of the Hamiltonian. These wavefunctionals 
are linear combinations of our basis operators: 



n n. 



®n = l^ b^Tm^) = ^ VkJ&nfa) (26) 

i i 

where the coefficients are chosen so that we have the normalisation condition 

(T4(0,(O))T4(^(O))) = 1 (27) 

The purpose of this common normalisation is to make a comparison of the coefficients b\ { for 
different 1Z, meaningful. 

The basic idea now is that the coefficients b\ { encode the 'shape' of the state corre- 
sponding to the wavefunctional, because they multiply the same operators, albeit in different 
representations, and with a common normalisation. So to compare our excited k = 1 and 
k = 2A wavefunctionals we need simply compare their coefficients and b 2Ai . This would 
be straightforward if the fa basis operators were orthogonal, but since they are not, and some 
have very large overlaps with each other, apparently different sets of b^s can in fact corre- 
spond to very similar wavefunctionals. Normally the obvious way to surmount such obstacle 
would be to calculate directly the overlap of the two wavefunctionals we want to compare. 
Unfortunately here that does not work because the operators have different iV-ality so that 
the overlap will be zero. So to proceed we need to transform the iV-ality of the k = 2A 
wavefunctional to k — 1 while preserving its essential spatial characteristics. We do so by the 
simple substitution (with t — throughout) 

^A = f^b n 2Al Tr' 2A (fa) — > = X^Tfyfc) ( 28 ) 

i i 

where we replace Tr' 2A by Tr'^ in the expression for ^ 2 a- Intuitively it is plausible to think of 
$2 A as having approximately the same 'shape' as ^ 2 a-> so ^ na ^ if we now compare $ 2A directly 
with any by the calculation of their normalised overlap 



< $"'t / $p' >V2< $2A%A > V2 



On,,n- w 1/0 . 2 ^„ t ~ . 1/0 (29) 



we are in fact comparing the shape of with that of our original & 2A wavefunctional. 

The above method of comparison clearly has a strong heuristic component and to some 
extent our confidence in the result will depend on how clear-cut it proves to be. So we take our 
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four lightest k = 2 A (variational) eigenstates, form the corresponding set of {3^4; ft = 0, 3}, 
and take the overlap of each one of these with the corresponding k = 1 (variational) eigenstates, 
{QJ;n = 0, ...,3}. Note that the former are not necessarily orthogonal, while the latter are, 
by construction (just as were the original {$2a})- We begin by showing the results, in Fig. [T5l 
for our longest / = 32a flux tube in SU(5). We see that the k = 2 A ground state has an 
overlap that is virtually 100% on the k = 1 ground state. Since we have already established 
that both ground states are to a very good approximation described by the free string theory, 
this result serves to give us confidence in this method of comparison. Turning now to the first 
excited k = 2A state, which is the one whose identity we are most interested in here, we see 
from Fig. [15] that this has an overlap that is almost entirely on the first excited k = 1 state. 
Since we have established that the latter is a Nambu-Goto-like stringy excitation, we infer 
that so is the first excited k = 2A state. Thus we have our answer: this state is tending to 
Nambu-Goto and is not an 'extra' state involving massive excitations of the flux tube. 

We also see from Fig. [15] that the next two P — +, q = k = 2A excited states are very 
similar to the corresponding k = 1 states and therefore also Nambu-Goto-like. In Fig. [TBI we 
show the corresponding plots for length I = 16a. The results for the ground and first excited 
states are unambiguous, but the comparison for the higher excited states is now less clear-cut. 
Perhaps not surprisingly, the states become more clearly Nambu-Goto-like as we go to longer 
strings. 

We have repeated this whole exercise for the k = 2S representation. Once again we find 
that the ground state has an almost 100% overlap onto the k = 1 ground state. However 
the first excited state is less clear-cut than for k = 2A, and the comparison for other states 
produces no useful information at all. 

We have learned that the first excited k = 2A state in the P = +, q = sector is in fact 
Nambu-Goto-like, so it has two phonons of momentum p = ±2tt/1. The fact that there is 
such a large correction to the free string energy when I is not very large, tells us that these 
phonons have a large interaction energy. On the other hand we have seen that the ground 
state k = 2 A state with P = +, q = 2 has an energy very close to the free string energy at 
all I. This state has two phonons of the same momentum p = 2-k/I. Thus we infer that the 
interaction energy of two phonons with the same momentum is very small. So although we 
have not found an 'extra' state, we have learned something very specific and interesting about 
the inter-phonon interactions. 

8 Conclusions 

In this paper we have examined to what extent k = 2 flux tubes that are wound around 
a spatial torus can be described by an effective string theory. This study complements our 
earlier study of fundamental {k = 1) flux loops [2], where we found that a simple Nambu-Goto 
free string description described the spectrum very well, even when the length of the flux tube 
was on the order of its width. Since one can regard a k = 2 flux loop as a bound state of 
two k = 1 flux loops, one would naively expect the extra internal structure to lead to larger 
corrections to any asymptotic free string description, and indeed that there should be extra 
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states that do not look like excited string states and which arise through the excitation of the 
massive binding modes. 

Because k = 2 states are much more massive than k — 1 states, some systematic errors 
are much more important for the former. This has meant that our most reliable calculations 
are for the ground states in the various sectors labeled by different quantum numbers. These 
quantum numbers are the parity P and the momentum along the flux tube, q. Since we expect 
the ground state unexcited flux loop to be invariant under translations along its axis, and to 
be symmetric under reflections across its axis, having P ^ + and/or q ^ explicitly probes 
the non-trivial excitation spectrum of the flux loop (in a way that transverse momenta do 
not, which is why we do not consider them). 

In addition to these exact quantum numbers, we have confirmed by explicit calculation 
that there is very little overlap between the totally symmetric (k = 2S) and totally antisym- 
metric (k = 2A) representations. We thus find that our variational eigenfunctionals fall into 
these representations to a good approximation, so that we can categorise the spectra by their 
behaviour under the full SU(iV) group, and not just its center Z^. 

In the k = 2 A sector we found that the ground states for q = 0,1,2 have energies that are 
very close to those of the free bosonic string theory (Nambu-Goto in flat D = 2 + 1). The 
leading permitted correction to (E(q = 0; 1) / (al)) 2 is 0(l/(y/al) 5 ), and its fitted coefficient 
has a natural value of 0(1). This is in contrast to the fundamental flux loop where the 
corresponding coefficient is unnaturally tiny, < 0(1/10). More surprising is that the q — 1,2 
ground state energies are well described by Nambu-Goto all the way down to l^fa ~ 1.5. Such 
a flux tube is a blob rather than a string and it is remarkable that the thin string oscillation 
spectrum should accurately describe its lowest modes. Just as for the k = 1 flux tube, one 
is tempted to read into this evidence for some kind of gauge-string duality that makes a thin 
string the starting point for describing even a short blob-like flux tube. For the low-lying 
excited states in the q = sector the deviations are much larger - and of a size that one would 
expect within an effective string theory description that was only valid asymptotically. 

The much heavier k = 2S sector is clearly afflicted by much larger systematic errors, which 
strongly limits our ability to draw any useful conclusions, although here too the ground states 
of various momenta q are recognisably string like. 

As an interesting aside, by including operators that wound twice around the torus, but 
such that the second winding overlapped little in space with the first, we were able to identify 
some eigenstates that correspond to doubly wound fundamental flux loops, in addition to the 
bound k — 2 singly wound flux loops. 

While it is very informative to observe how well the various states are described by simple 
free string excitations, this can only be the first step. We also need to pin-point the nature 
of the corrections to the free string theory. This involves not only quantifying and searching 
for patterns in the corrections to the free string excitations - the interactions amongst the 
phonons - but also identifying excited states that are additional to those associated with 
simple string oscillations. Such additional states, due to excitations of massive modes, would 
normally be predicted by any theoretical approach. 

For example, in a naive field theoretic approach, where the flux tubes are some kind of dy- 
namically generated non-Abelian Nielsen-Olesen vortices, there should be massive excitations 



21 



associated with the nontrivial flux tube structure. In an 'AdS/CFT' approach there would be 
a highly curved metric in the region of the higher dimensional space where the string picks up 
its non-trivial energy per unit length. Naively, the energy of such an 'extra' state should be 
some (roughly) constant amount /i above the ground state energy. As we remarked earlier, a 
naive estimate for // might be the lightest glueball mass, ~ 4-^0"/, or possibly the gap between 
the lightest glueball masses, ~ On such an estimate, such states should be easily visible 

for small flux tube lengths, where they would be lighter than some of the simple thin string 
excitations. A glance at the q = spectra in Figs. [5] and El immediately suggests that the first 
excited state could be a candidate for such a state. However our analysis in Section [7] of the 
wavefunctional of that state showed that its 'shape' was too much like that of the correspond- 
ing k = 1 flux tube, which is unambiguously Nambu-Goto-like, to believe that it is anything 
other than a simple string excited state. We have in fact not observed any states, in the 
range of energies where we have some control, that could be interpreted as such non-stringy 
excitations. This could either mean that we have a poor overlap onto such states or that their 
masses are much larger than expected. (The glueball mass might not be a good guide if its 
dynamics is that of a string-like contractible closed loop of flux, so that neither its mass nor 
the gap to the first excited state are strongly influenced by the dynamics associated with the 
non-stringy massive modes.) In the former case it means that we must suitably extend our 
basis of operators; and for that purpose it would be useful to have some guidance from, say, 
the gauge-gravity side, about the lightest such extra states. In the latter case, the implication 
would be that the non-stringy physics was essentially decoupled from the low-energy physics 
of the confining gauge theory. Resolving this issue would provide something very interesting 
in either case. 

Although the question of non-stringy 'extra' states has not been settled in this paper, we 
have been able to identify some qualitative features of the interactions amongst the phonons 
along the flux tube. As we have just remarked, we have shown that the first excited q = 
P = + state is in fact a string excited state, i.e. it has two phonons of opposite momentum 
p = ±2tt/1. And, as we see in Figs. [5] and El the energy of this state has very large corrections 
to the corresponding free string energy. On the other hand the lightest q = 2 P = + state, 
which has exactly the same two phonons, but with the same momentum, has as we see in 
Figs, d] and [21 an energy that is virtually identical to that of the corresponding free string 
state, over our whole range of I. Thus we can come to a very specific conclusion about the non- 
trivial interactions between two phonons: they are negligible when the phonons have exactly 
the same momenta (and are therefore at 'threshold'), but are large when the momenta are 
exactly opposite. 

One might ask why such a conclusion had not already been reached in our much more 
accurate calculations of the fundamental k = 1 flux loop spectrum in [2]. There are at least 
two plausible reasons. Firstly, we have seen that the corrections are in general much smaller 
for k = 1 than for k = 2 flux loops, so much so that the first excited q = 0P = + k= l flux 
loop looks, at a first glance, to be close to the free string prediction. Secondly, we observed 
in [2j some significant corrections in the q = 2 spectrum at small I, in contrast to the case 
of k = 2. However we have since been able to show that these are largely due to lattice 
spacing corrections to the energy-momentum dispersion relation. Thus a renewed analysis 
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may well show that the k — 1 flux loop displays the same kind of non-trivial interaction 
between phonons that we have seen with the k = 2 loop. The fact that the k = 2 calculation 
in this paper has been performed at a smaller value of a, together with the fact that the 
k = 2 loops are much heavier and therefore less sensitive to the momentum and to any lattice 
corrections thereof, may be part of the reason that it has been easier to identify this effect in 
the present, nominally more difficult, calculation. 

There is of course much more information about the phonon interactions, that is implicit in 
our calculations. We have focused on the above single result because it involves energy levels 
that are non-degenerate, making the argument and conclusion particularly straightforward. In 
most other cases corrections to the Nambu-Goto energy levels are accompanied by splittings 
and (presumably) mixing of the would-be degenerate string states, and this will complicate 
the analysis. While there are other examples of non-degenerate states that can be usefully 
compared, e.g. the ground state in the P = +, q = 3 sector and the first excited state in the 
P — +, q — 1 sector, these involves values of momentum q > 2 which have not been studied 
in the k = 2 calculation of the present paper (but have been included in our accompanying 
k = 1 calculations). We therefore leave a more complete analysis to a forthcoming publication 
[T7] that will describe in detail all our k = 1 and k = 2 calculations. 
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Appendix A: Construction of operators 



To calculate the excitation spectrum of k = 2 strings, it is necessary to find a way to project 
onto such states. The way to achieve this is to find a suitable basis described by good quantum 
numbers, in which our operators will be encoded. In our case, this basis is defined by the 
quantum numbers of momentum, q, along the string, and the parity P, where the latter tells 
us how a string transforms under reflections in its transverse axis. This imples that we need to 
introduce transverse deformations in Polyakov loops and construct line paths that transform 
in a certain way under such parity reflections and longtitudinal momentum. In general, the 
more operators we use the better the results we obtain; we therefore construct a plethora of 
Polyakov paths trying to extract states with high overlaps. 

In the first step of the calculation we construct the most commonly used k = 2 string 
operators: (Tr{l p }) 2 and Tr{/p}. Our desire to extract the excitation spectrum and introduce 
more degrees of freedom, such as the transverse parity and momentum, forces us to modify 
them and thus to complicate their structure, as discussed in Section 12.31 In Eq. (13011321) we 
demonstrate how these new operators are composed for the k = 1, 2 cases. For simplicity, be- 
low we demonstrate the way our operators have been constructed with a particular transverse 
deformation that is easy to visualize. First are the negative and positive parity operators for 
the k = 1 strings: 

C? = Tr { _TL_ } ± ~LT- } • (30) 

From these we construct the two simplest sets of positive and negative parity operators for 
the k = 2 string, i.e. those corresponding to Tr{l p } 2 

$Lt = Tr{ ^} 2 ±Tr{^} 2 , (31) 

and those corresponding to Tr {l 2 } 

<Ct=Tr { _TL_ ■ _TL_ } ± ~LT- ~LT- } • (32) 

Next, we project onto the k = 2 totally antisymmetric and symmetric representations. 
To single-out the irreducible representations which describe the theory, we perform antisym- 
metrisation and symmetrisation according to the relevant Young-tableau decomposition of 
k = 2 fundamental colour sources: □(££)□= I ® H- The resulting operators will be of 
the type: Tr{/ P } 2 ± Trj/^}, where +(— ) for symmetric (antisymmetric) representation. Once 
more we need to introduce transverse deformations, in order to project onto the non-trivial 
irreducible representations that characterise the closed flux tube in D = 2 + 1. Examples of 
such operators are demonstrated in Eqs. (1331 1341) below. Begining with the projection onto 
the k = 2 antisymmetric representation we find 

$£=± = [Tr { _TL_ } 2 - Tr { _TL_ ■ _TL_ }] ± [Tr { -"LT f - Tr { ~u~ ■ ~VT }], (33) 
and projecting onto the k = 2 symmetric representation one obtains 
C = 25 = [Tr {^} 2 + Tr { _TL_ ■ _TL_ }] ± [Tr { ""LT f + Tr { -LT" ■ ~U~ }]. (34) 
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As the notation suggests, in both cases above the ± signs determine the parity P — ± of 
The complete set of polyakov lines used in our calculation is presented in Table HI 

As mentioned in Section 12.31 we have also attempted to test whether the k = 2 string 
spectrum includes unbound w = 2 states, which we expect will only appear in our calculation if 
we use the appropriate operators. These new states are expected to be described by frequencies 
lower than those describing the w = 2 bound states. Although some k = 2 operators look as 
if they wind twice around the torus i.e Eqs. (I3TH32]) . the way we had been constructing them 
prohibits us to project onto states with lower frequencies since each Polyakov loop starts and 
ends at the same lattice point within one lattice size. To overcome this we construct Polyakov 
lines that wind twice around the torus with transverse deformations at the joint of the two 
lattices as in : 

C= 2 ± = Tr { n }±Tr{ U }. (35) 

Note the number of lattice links in the x directions that goes into each of the operators in the 
traces above is 2L X , i.e. twice the lattice size. This is in contrast to the case of Eqs. fl3T|) - fl32|) . 
where it is only L x . 

Finally, let us note that the new operators in Eq. [3S1 transform the same under the centre 
of the group Z N as the other k = 2 operators i.e Q w =2 = Tr(/p/' p ) — > z 2 Tr(l p l' p ). So these new 
'unbound' operators have A-ality k = 2 and will therefore contribute to an extended k = 2 
spectrum. 
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Table 1: The seven lowest Nambu-Goto energy levels for the w = 1 closed string. If the 
number of creation operators is even(odd) then the state has positive (negative) parity. 
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Table 2: Overlaps of Polyakov loops in the k=2A and k=2S representations, with blocking 
levels b A and b s respectively, at t = and as defined by eqn (!25l) . For SU(4) at /3 — 32.0 on 
I x L x 20 lattices. 
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Table 3: Overlaps of Polyakov loops in the k=2A and k=2S representations, with blocking 
levels b A and bs respectively, at t — 0, 1, 2 and as defined by eqn ([2"5"j) . For SU(5) at (3 = 80.0 
on 24 2 32 lattices. 
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Table 4: The sixteen transverse deformations used in the construction of the operators in 
this work. Each line comes in five different blocking/ smearing levels, in order to enhance the 
overlap onto the physical states. 
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Figure 1: Energies of the lightest four states with q — 0,1, 2, for winding flux tubes in the k = 2 
antisymmetric representation of SU(4) at (5 = 50. Lines are the Nambu-Goto predictions of 
eqn (IT5|) . with the string tension obtained by fitting the ground state. 
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Figure 2: Energies of the lightest four states with q = 0,1,2, for winding flux tubes in the k = 2 
antisymmetric representation of SU(5) at (3 = 80. Lines are the Nambu-Goto predictions of 
eqn (fn?|) . with the string tension obtained by fitting the ground state. 
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Figure 3: Flux tube excitation energies from the SU(4) q = 1,2 calculations in Fig.[TJ extracted 
using eqn (!2ip . Lines are the Nambu-Goto predictions from eqn (fT3|) . 
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Figure 4: Flux tube excitation energies from the SU(5) q — 1,2 calculations in Fig. [21 extracted 
using eqn( j2il . Lines are the Nambu-Goto predictions from eqn(fl~3l. 
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Figure 5: Energies of the lightest four states with P = + and g = 0, for winding flux tubes 
in the k = 2 antisymmetric representation of SU (4) at (3 = 50. Lines are the Nambu-Goto 
predictions of eqn (fT3|) . with the string tension obtained by fitting the ground state. 
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Figure 6: Energies of the lightest four states with P = + and q = 0, for winding flux tubes 
in the k = 2 antisymmetric representation of SU(5) at /3 = 80. Lines are the Nambu-Goto 
predictions of eqn ffl3|) . with the string tension obtained by fitting the ground state. 
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Figure 7: Energies of the lightest four states with q = 0, 1,2, for winding flux tubes in the 
k = 2 symmetric representation of S'6 r (4) at /? = 50. Lines are the Nambu-Goto predictions 
of eqn (fT3|) , with the string tension obtained by fitting the ground state. 
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Figure 8: Energies of the lightest four states with q = 0, 1,2, for winding flux tubes in the 
k = 2 symmetric representation of SU(5) at /3 = 80. Lines are the Nambu-Goto predictions 
of eqn( ri"3l . with the string tension obtained by fitting the ground state. 
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Figure 9: Flux tube excitation energies from the SU(5) q = 1, 2 calculations in Fig. [HI extracted 
using eqn fpZTI) . Lines are the Nambu-Goto predictions from eqn (fl~3"]) . 
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Figure 10: The overlap \OAs(hj',t = 0)1 as defined by eqn fl25l) for SU(A), L = 32a, 
and q = with i, j = 1 — 80. 
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Figure 11: Overlaps squared onto k = 2 A (red circle) and k = 2S (blue square) of the low- lying 
k = 2 flux tube states with the lengths I shown. In sectors, from left to right, {P — +, q — 0}, 
{P = +, q = 0}, {P = -, q = 0}, {P = -, q = 1}. All in SU(4) at (3 = 50. 
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Figure 12: Spectrum of = 2 flux loops with P = — and g = 0, as extracted using the 
extended basis of operators, that includes w = 2 k = 1 operators as in eqnfjll). 
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Figure 13: Spectrum of k = 2 flux loops with P = — and g = 0, as extracted using only our 
usual normal k = 2 basis of {£; = 2A} @{&; = 2S 1 } operators, as in eqn([3]), and excluding the 

oj = 2 states 




Figure 14: The contribution of each of our extended basis of operators to the ground and first 
excited states shown in Fig. [12] for L = 16a. For i < 100 the operators are the ordinary k = 2 
ones and for i > 100 the operators are the additional w = 2 ones. 
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Figure 15: The Overlaps defined in Eq. (l29l) for the lightest four states in the SU(5) k 
antisymmetric representation with P = +, q = and flux tube length I = 32a. 
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Figure 16: The Overlaps defined in Eq.( l29i) for the lightest four states in the SU(5) k 
antisymmetric representation with P = +, q = and flux tube length I = 16a. 
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